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We show the possibility of extracting important information on the symmetry term of 
the Equation of State (EOS) directly from multifragmentation reactions using stable iso- 
topes with different charge asymmetries. We study n-rich and n-poor Sn + Sn collisions 
at 50AMeV using a new stochastic transport approach with all isospin effects suitably 
accounted for. For central collisions a chemical component in the spinodal instabilities is 
clearly seen. This effect is reduced in the neck fragmentation observed for semiperipheral 
collisions, pointing to a different nature of the instability. In spite of the low asymmetry 
tested with stable isotopes the results are showing an interesting and promising depen- 
dence on the stiffness of the symmetry term, with an indication towards an increase of 
the repulsion above normal density. 

1. Introduction 

Our starting point is that the key question in the physics of unstable nuclei is the knowl- 
edge of the EOS for asymmetric nuclear matter away from normal conditions. We remark 
the effect of the symmetry term at low densities on the neutron skin structure, while the 
knowledge in high densities region is crucial for supernovae dynamics and neutron star 
cooling [0-0. 

Effective interactions are obviously tuned to symmetry properties around normal con- 
ditions and any extrapolation can be quite dangerous. Microscopic approaches based on 
realistic NN interactions, Brueckner or variational schemes, or on effective field theories 
show a quite large variety of predictions, see Fig.l. 

In the reaction dynamics with intermediate energy radioactive beams we can probe 
highly asymmetric nuclear matter in compressed as well as dilute phases: the aim of this 
paper is to show that fragmentation events have new features due to isospin effects and 
that some observables are particularly sensitive to the symmetry term of the EOS. 

For fragmentation mechanisms a new qualitative feature is expected: the onset of chem- 
ical effects on spinodal instabilities [jSHi~0|l, due to the intrinsic coupling between isoscalar 
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Fig.l - EOS for various effective forces, SKM* [§ (dashed), SLy230b § (squares) and 
BPAL32 pF (crosses). Top: neutron matter (up), symmetric matter (down); Bottom: po- 
tential symmetry term. 



and isovector modes. The formation of a very neutron rich gas vs. an almost symmet- 
ric liquid phase is expected for neutron excess systems, in a dynamical non-equilibrium 



mechanism on short time scales [10-12]. This effect, which appears very much reduced 
in a statistical multifragmentation calculation |Tj|, will show up in a production of more 
stable primary intermediate mass fragments and in enhanced yields for neutron rich light 
isobars. 

All these predictions were actually based on linear response approaches [|10],[lTJ, nuclear 



matter dynamics in a box [|TT] , |T2|[ or on the evolution of a suitably prepared excited source 
131. In this paper we will present the first results of a fully "ab initio" calculation of 



fragmentation reactions for systems with different charge asymmetry. 



2. "Ab initio" simulations of Sn + Sn reactions 

A new code for the solution of microscopic transport equations, the Stochastic Iso-BNV, 
has been written where asymmetry effects are suitably accounted for and the dynamics of 



fluctuations is included fl4| , [T5|j . A density dependent symmetry term is used also in the 
ground state construction of the initial conditions. Isospin effects on nucleon cross section 
[|T6f and Pauli blocking are consistently evaluated. In order to simplify the analysis of the 
most sensitive observables to isospin effects we have chosen a Skyrme force with the same 
soft EOS for symmetric Nuclear Matter (NM) (K = 201MeV) and with two different 
choices for the density dependence of the symmetry term, see Fig.l (bottom), asy-stiff 
(like in BPAL32, see also [17]) and asy-soft (like in SKM*). In this way we force the 
symmetric part of the EOS to be exactly the same in order to disentangle dynamical 
symmetry term effects. We will show that the reaction mechanism is sensitive to the 



different behaviours, although with stable nuclei we will limit the possible asymmetries 
and moreover we will certainly not reach high compression regions. 

We have studied the 50AMeV collisions of the systems 124 Sn+ 12A Sn and ll2 Sn+ 112 Sn, 
where new data are under analysis at NSCL — MSU ||19| . We will comment 100 events 
generated in semi-central (b = 2fm) and semi-peripheral (b = 6/m) reactions. In Fig. 2 we 
show a typical impact parameter evolution of the density plot (projected on the reaction 
plane) for one event (neutron rich case, asy-stiff EOS). 
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Fig. 2 - l2A Sn -\- 12i Sn 50AM eV: time evolution of the nucleon density projected on 
the reaction plane. First two columns: b = 2fm collision, approaching, compression and 
separation phases. Third column: b = Afm , fourth column: b = 6/m, separation phase up 
to the freeze- out. 



We remark: i) In the cluster formation we see a quite clear transition from bulk |2C 22 

ii) The 



to neck |23|-^| instabilities, ii) The "Freeze- Out times", when the nuclear interaction 
among clusters disappears, are decreasing with impact parameter. These two dynamical 
effects will influence the isospin content of the produced primary fragments, as shown 
later. 

A detailed analysis of the results from 100 events for the same system {asy-stiff calcu- 
lation) is shown in Fig. 3 (central, b = 2/m) and Fig.4 (peripheral, b = 6/m). Each figure 
is organized in this way: 

Top row, time evolution of: (a) Mass in the liquid (up) and gas (down) phase; (b) 
Asymmetry I = (N — Z) / (N + Z) in the gas "central" (solid line and squares), gas total 
(dashed+squares), liquid "central" (solid+circles) and clusters (stars). "Central" means 
in a cubic box of side 20/m around the cm.. The horizontal line shows the initial average 
asymmetry; (c) Mean Fragment Multiplicity Z > 3. The saturation of this curve defines 
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Fig. 4 - Like Fig.3 for b = 6/m. 



the freeze-out configuration, as we can also check from the density plots like in Fig. 2. 

Bottom row, properties of the "primary" fragments in the Freeze- Out Configuration: 
(d) Charge Distribution, (e) Asymmetry Distribution and (f) Fragment Multiplicity Dis- 
tribution (normalized to 1). 
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Fig. 6 - Like Fig.5 for b — 6/to. 



We see a neutron dominated prompt particle emission and a second neutron burst at 
the time of fragment formation in the "central region". The latter is consistent with the 
dynamical spinodal mechanism in dilute asymmetric nuclear matter, as discussed before. 
The effect is quite reduced for semi-peripheral collisions (compare Fig.s 3b and 4b) and 
the IMF's (3 < Z < 12) produced in the neck are more neutron rich (Fig.s 3e and 4e). 
This seems to indicate a different nature of the fragmentation mechanism in central and 
neck regions, i.e. a transition from volume to shape instabilities with different isospin 
dynamics. In more peripheral collisions the interaction time scale is also very reduced 
(Fig. 4c) and this will quench the isospin migration. 

In Fig.s 5,6 we have the corresponding results for the asy-soft EOS. The main qual- 
itative difference is a larger prompt neutron emission (compare the gas asymmetry in 
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Fig. 7 - Correlation between mean IMF multiplicity and charge of the heaviest fragment: 
squares n-rich " asy-stiff" ; diamonds n-poor " asy-stifF' ; circles n-rich " asy-soft" ; stars n-poor 
" asy-soft" . 



Fig.s 3b,4b with the corresponding Fig.s 5b, 6b) joined to a "slower" dynamics (see the 
(c)-plots), as expected from the more attractive nature of the asymmetric EOS fL2f . 

The interesting point is that the effect on fragment production is different for central 
and peripheral events. With respect to the asy-stiff case we have a smaller mean IMF 
multiplicity for central collisions (Fig.s 3f vs. 5f) and more fragments produced in the 
neck region (Fig.s 4f vs. 6f). 



Neck instabilities have also the feature of forming clusters on the "spectator" side |23 
leading to a "fission-like" splitting of the Projectile-/Target-like fragment. This mecha- 
nism has been clearly observed in accurate kinematical selections [p6fl . In our study of 



n-rich systems the rate of such dynamical fission processes is systematically larger for 
the asy-soft EOS, as expected from the previous discussion. This seems to be an quite 
sensitive observable to look at. 

The most important qualitative difference of the neutron poor case, U2 Sn + 112 Sn, 
is a larger prompt proton emission, in particular for central collisions. The number of 
protons available to produce clusters is then reduced with a related smaller mean IMF 
multiplicity. Moreover the reduction of the overall asymmetry is acting in the same 
direction. In Fig. 7 we show the average IMF multiplicity vs. the average charge of the 
heaviest produced fragment, which measures the centrality of the collision, for the two 
symmetry terms and for the two systems. The effect, in agreement with recent data pTf , is 
more evident with the asy-stiff choice, in particular for central collisions (low < Z >heavy)- 

3. Outlook 

Starting from simulations of reaction dynamics performed with a new transport code 
where isospin and fluctuation effects are suitably accounted for, we have shown that dis- 
sipative heavy ion collisions, and in particular fragment production reactions, at medium 
energies are rather sensitive to the density dependence of the symmetry contribution to 
the nuclear Equation of State. 

Interesting differences appear between isospin effects on central and neck fragmentation. 



An indication for a "stiff' symmetry term is emerging, consistent with transverse flow 



calculations in the same energy range 28 



The effects are not large but quite encouraging for similar studies with radioactive 
beams. A possibility is emerging of obtaining in terrestrial accelerator laboratories im- 
portant information on the symmetry term of large astrophysical interest. It appears 
essential to have good charge asymmetric beams available at intermediate energies and 
to perform more exclusive experiments. 
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